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Abstract 

The generalized tight-binding model is developed to investigate the feature-rich 
magneto-optical properties of AAB-stacked trilayer graphene. Three intragroup and 
six intergroup inter-Landau-level (inter-LL) optical excitations largely enrich the 
magneto-absorption peaks. In general, the former are much higher than the latter, 
depending on the phases and amplitudes of LL wavefunctions. The absorption spec¬ 
tra exhibit the single- or twin-peak structures which are determined by the quantum 
modes, LL energy spectra and Fermion distribution. The splitting LLs, with different 
localization centers (2/6 and 4/6 positions in a unit cell), can generate very distinct 
absorption spectra. There exist extra single peaks because of LL anti-crossings. AAB, 
AAA, ABA, and ABC stackings quite differ from one another in terms of the inter-LL 
category, frequency, intensity, and structure of absorption peaks. The main charac¬ 
teristics of LL wavefunctions and energy spectra and the Fermi-Dirac function are 
responsible for the configuration-enriched magneto-optical spectra. 


Keywords: AAB-stacked trilayer graphene, magneto-absorption spectra, optical se¬ 
lection rules 


* Corresponding author. Tel: -1-886-6-275-7575. 

E-mail addresses: 128981084@mail.ncku.edu.tw (C.Y. Lin), mflin@mail.ncku.edu.tw 
(M.F. Lin) 


1 



1. Introduction 


Few- and multi-layer graphenes with different stackings have attracted a lot of experi¬ 


mental lH9 and theoretical 10 -20 studies because of the special hexagonal symmetry. 


They have been successfully synthesized by experimental methods, such as exfoliation 
of highly orientated pyrolytic graphite 0i , chemical and electrochemical reduction of 
graphene oxide [^|^, metalorganic chemical vapour deposition (MOCVD) [5]-[^, and arc 
discharge [^. Graphene shows unusual essential physical properties, including the di¬ 
verse optical selection rules 10 rich magnetic quantization 10 15 , half-integer Hall 


effect [3,16-18 , low frequency plasma 19,20 , and others. Such properties could be easily 


modulated by stacking configurations 14| 21|-|23 , dopping 24| 25 , layer number 26| 27 
magnetic field (B = B^z) 


29 


, electric field 


30 


32 


, mechanical strain 


33 


35 


, and 


temperature 36,37 . To date, four important typical stacking configurations, AAB 


ABA [4]-[^, AAA and ABC [^|^, have been discovered. In this work, we develop 

the generalized tight-binding model to investigate the magneto-absorption spectra of AAB- 
stacked trilayer graphene. A detailed comparison among four stacking trilayer systems is 
also made. 

The rich magnetic quantization is induced by the distinct stacking configurations. A 
N-layer graphene possesses N groups of LLs, where each group consists of valence and 

Monolayer graphene has the A/Ho-dependent inter-LL transition 


conduction ones 


14 


energies, in which the absorption peaks satisfy the selection rule of An = ±1 43 . Several 


theoretical studies on the magneto-optical properties have been elaborated for the AAA- 
and ABC-stacked trilayer graphenes. The AAA stacking presents three pairs of Dirac 


cones results in only three intragroup inter-LL transition categories 38,39 . The above 
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mentioned \fBo dependence and selection rules are also valid for N-layer AAA graphene. 
The ABC stacking has three pairs of energy bands, including the linear, sombrero-shaped. 


and parabolic bands, which give rise to nine categories of inter-LL excitations 40,41 . On 


the other hand, a thorough systematic study on the magneto-optical properties of the ABA- 


stacked trilayer graphene is absent up to now 41,42 . The electronic structure exhibits two 


pairs of parabolic bands and a pair of linear ones. This means that the ABA stacking could 
be regarded as the combination of monolayer and AB-stacked bilayer graphenes in optical 
properties. It is worth mentioning that the AB bilayer system presents four categories of 
absorption peaks, in which there is no simple relationship between the inter-LL transition 


energies and Bq 43,44. The ABA graphene is thus expected to have hve categories of 


inter-LL excitations, what will be investigated and presented in this work (Figs. 6(a)-6(c)). 
The energy dispersions of AAB-stacked trilayer graphene consist of three divergent pairs 
of valence and conduction bands, oscillatory (5’^’’^), sombrero-shaped and parabolic 

ones ((Fg(Fig. 1(b)) [^. Signihcantly, the bands (Fig. 1(c)) are hrst revealed in the 
AAB stacking implying that its electronic structure is exclusive among the aforementioned 
trilayer systems. How the unique electronic properties reflected in the magneto-absorption 
spectra is worthy of a thorough study. 

The methods of optical spectroscopies are powerful tools to identify the optical exci¬ 
tations and selection rules in condensed-matter systems. Such experimental techniques, 
including absorption, transmission, reflection, Raman scattering and Rayleigh scattering 
spectroscopies, have been successfully applied to two-dimensional (2D) materials, concern¬ 


ing a large class of graphene-related systems, such as, few-layer graphenes 45 , carbon nan¬ 


otubes 46 , graphene nanoribbons 47 , graphite and graphite intercalation compound 48 ; 
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also other few-layer systems, e.g., silicenes 49,50 , topological insulators 51 , M 0 S 2 52,53 , 


and so on. The intra- and inter-LL optical transitions in graphene systems have been ob¬ 
served in infrared spectroscopy. Evidently, the held-dependent frequency and the selection 


rule of An = ±1 are identihed by transmission and Raman scattering spectroscopy 54 -59 . 


In particular, the square-root and linear Ro-dependent frequencies due to the massless and 
massive Dirac-fermions, respectively, of the monolayer and AB-stacked bilayer graphenes. 


are revealed in the infrared transmission spectra 54 56 . These Dirac-fermions are further 


conhrmed by the magneto-Raman spectroscopy for the few-layer AB-stacked and ABC- 


stacked graphenes 57-59 . It should be noted that the experimental measurements men¬ 


tioned previously are only resolved on the intragroup magneto-absorption spectra due to 
the hrst group of LLs. 

The center of interest in this work is to investigate the conhguration-enriched magneto¬ 
optical properties of AAB-stacked trilayer graphene using the generalized tight-binding 
model. The special electronic properties are directly reflected on the LL energies and 
wavefunctions, so that they induce the extraordinary features in the magneto-optical prop¬ 
erties. This study shows that there exist nine categories of absorption peaks, including three 
intragroup and six intergroup ones. The intragroup transitions remarkably differ from the 
intergroup ones in the intensity, energy spacing, and structure of peaks. The former ex¬ 
hibits very low intensity, while the latter are much higher and comparable with those of 
monolayer graphene. Moreover, the splitting LLs, arising from the destruction of inversion 
symmetry, can diversify the magneto-absorption spectra. These characteristics are absent 
in the other trilayer graphene systems. There are important differences among four typical 


stacking trilayer graphenes in the excitation category and the main features of absorption 
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peaks. The predicted results, with regard to nine categories of inter-LL transitions, can be 


verihed by optical spectroscopies 54 59 
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2. Method 


The generalized tight-binding model is developed to study the magneto-electronic spec¬ 
tra of AAB-stacked trilayer graphene [^. For the hrst two layers, shown in Fig. 1(a), 
all carbon atoms have the same (x,y) projections while the third layer can be obtained 
by shifting the hrst (or the second) one by a distance of b along the armchair direction. 
The parameters of ten intra- and inter-layer atomic interactions are used (70 = —2.569 
eV, 7 i = -0.263 eV, 72 = 0.32 eV, 73 = -0.413 eV, 74 = -0.177 eV, 75 = -0.319 eV, 
76 = —0.013 eV, 77 = —0.0177 eV, 73 = —0.0319 eV, and 79 = —0.012 eV ) to characterize 
the special energy bands. A uniform perpendicular magnetic held can induce the Peierls 
phase related to the vector potential and thus enlarge the unit cell to become a rectangle 
including 127?^ carbon atoms {Rb = hc/eBo) (Fig. 1(a)). The eigenfunctions of the quite 
large Hamiltonian is solved by an exact diagonalization method to be in the arrangement 
of 12Rb tight-binding functions as follows. 


3 2Rb 


I'l-r) “EE 


l=l m=l 


Ca‘ 


m,k 




m,k 


+ Cb 


m.k 


tdI 

-°m,k 


in which c and v denote, respectively, the conduction and valence states; (cbi ) 

represent the subenvelope functions on the six distinct sublattices. For each LL, the quan¬ 
tum number (n^’*^) could be determined by the number of zero points in the dominating 
subenvelope wavefunctions. 

When ABA-stacked trilayer graphene is present in an electromagnetic held with the 


electric polarization E \\ y, a.n occupied LL will be excited to an unoccupied one. The 
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optical absorption function, according to the Fermi golden rule, is expressed as 


A{uj) oc / 
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-£’'=(k, m') — E’'(k, m) — uj — iV y 
where P is the momentum operator, /(E'^’'^(k, m)) the Fermi-Dirac distribution function; 
F 1 meV) the broadening parameter. The absorption spectrum is associated with the 
velocity matrix elements (the hrst term) and the joint density of states (the second term). 
The former can determine whether the inter-LL transitions are available. The velocity 


matrix elements, as successfully done for carbon-related materials 60 , are evaluated under 


the gradient approximation in the form of 


T^(k,m') 
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Since the nearest-neighbor atomic interaction is much larger than the interlayer ones, only 
the Hamiltonian elements of the former dominate the velocity matrix elements. On the 
same layer, dH/dky is almost the same for all the nearest neighbors in the enlarged unit 
cell. Consequently, the magneto-optical selection rule depends on the relations between 
the initial state on the A^ {B^) sublattice and the final state on the (H^) one. That is, 
these two states on the different sublattices of the same layer must have the same quantized 


mode. 
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3. Magneto-absorption spectra 


AAB-stacked trilayer graphene exhibits three groups of LLs, where each group con¬ 
sists of occupied valence and unoccupied conduction ones. Each LL is characterized by 
the subenvelope functions on the six sub lattices, as indicated in Figs. 2 and 3 for the 
[kx = 0,/Cy = 0) state at 2/6 localization center under Bq = 40 T. The hrst, second and 
third groups are initiated at 0 (0), -0.23 (0.24) and -0.58 (0.55) eV for the valence (con¬ 
duction) states, as indicated by blue, green and red colors, respectively. The valence and 
conduction LLs are asymmetric about Ep = {] because of the interlayer atomic interactions. 
It is expected to induce the twin-peak structures in the magneto-absorption spectra. The 
quantum numbers of the hrst ( second (n^’’^) and third groups {n^^) are, respectively, 
determined by the A^, and sublattices with the dominating subenvelop functions 
(Figs. 2 and 3). As to the hrst group, there exist three LLs near Ep, in which the middle 
= 0 LL is located right at Ep, while the other two are assigned = 1. The next 
LLs = 2, = 3, ... are in an usual monolayer-like ordering. Distinctly, the quantum 

numbers of begin from = 1, then = 0, = 2, = 3, and so on. As for 

it normally increases from zero for the higher conduction states, and the lower valence 
states. 

The stacking conhguration of AAB results in the special relationship between the diher- 
ent subenvelop functions. All the LL wavefunctions are almost well-behaved in the spatial 
distributions at Bq = 40 T, in which the main modes are much stronger than the side ones. 
The subenvelope functions of and {B^ and B^) are similar in the number of zero 


points since the hrst and second layers have the same (x,y) projection 15 . That is to say, 
the quantum numbers determined by the A^ and A^ sublattices are identical. It should be 



noted that, for the first and second groups (the third one), the phases of the and 
sublattices are identical for all the conduction (valence) states but opposite for the valence 
(conduction) ones. Moreover, the spatial distribution of certain subenvelope functions are 
much smaller compared to the others. For example, for the hrst group, the amplitude of 
the A^ sublattic is much weaker than those of the other sublattices (Fig. 2). The critical 
differences in amplitude and phase are expected to induce the unique absorption spectra. 

The nine categories of optical excitations, which include three intragroup and six inter¬ 
group ones, are clearly revealed in the magneto-absorption spectrum (Figs. 4(a)-4(e)). The 
inter-LL optical excitations between any two groups are available only under the require¬ 
ment of a specihc selection rule. The optical absorption peaks corresponding to the first, 
second, and third LL groups are respectively indicated by black, green, and red numbers 
(underline for valence LLs). The intragroup and intergroup absorption peaks quite differ 
from each other in the spectral intensities, peak structures, and energy spacings in twin 
peaks. The whole intragroup categories, excluding the threshold peak 01, are twin peaks 
with very low intensity. While for the intergroup ones, except for the transitions between 
the first and third LL groups, all of peaks show the single-peak structure with monolayer¬ 
like intensity. The intensity difference arises from the special relations of the subenvelope 
functions between the A’' sublattice of the initial states and the [A'') sublattice 
of the hnal ones (discussed later). Regarding the energy spacings in twin peaks, those of 
the intergroup excitations are much wider than the intragroup ones, mainly owing to the 
asymmetry of LL energy spectra. 

The intensities of absorption peaks are determined by the summation of {A^\B^) and 
(R^|W), where I = 1,2,3. As to the hrst and third groups, {A^\B^) ^ 0, especially for 
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low-lying LLs because of the very low amplitudes of the and subenvelope functions, 
respectively. Besides, the contributions of are fainter for higher or deeper LLs due 

to the lower amplitudes on two sublattices. Moreover, the peak intensity can be largely 
enhanced or reduced depdending on the phase relationship between the two sublattices: 
A^ {B^) of the initial states and B^ {A'') of the hnal ones. In general, the intensities of 
the intergroup transition peaks are much higher than those of the intragroup ones. For 
example, the intragroup peak 12 at a; ~ 0.06 eV (marked by blue arrow) is about 1/15 
compared with the intergroup one 1 1 at cu ~ 0.25 eV (marked by green arrow in Fig. 4(a)). 
For the former, the \B^) of the n\ = 1 initial state and the of the nl = 2 hnal state 
have the distinct phases, while \B^) and {A^\ possess the same ones. As a result, 
and {A^\B'^) present the opposite signs which reduce the velocity matrix elements, and 
thus lower the absorption peak intensity. On the contrary, the peak intensity of the latter 
is enhanced due to the fact that all the dominating subenvelope functions of the = 1 
initial and n\ = 1 hnal states on the sublattices, {B^, A^), {B‘^, A^); (A^, B^), have similar 
phase relations. It should be noted that there are a number of extremely low peaks in the 
spectra, e.g., 51, 51 (red arrow), and so on, due to the appearance of the side modes. Those 
transitions do not satisfy the optical selection rule and thus could be neglected. 

The intragroup inter-LL excitations possess a lot of low-intensity absorption peaks, as 
shown in Fig. 4(a). The threshold peak of n\n1 = 01, with frequency of a; = 0.01 eV, 
possesses an ignorable intensity (inserted in Fig. 4(a)). This is caused by the extremely 
weak subenvelope functions of = 0 LL in the B^, A^ and A^ sublattices. The other low 
peaks belong to the twin-peak structures, such as = 21 and = 12 (a; 0.05 

eV), n\nl = 32 and = 23 (a; ~ 0.13 eV), and so on. The quantum mode of the 
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n\ = 2 A^{A‘^) sublattice is identical to that of the = 1 one, so that the 21 peak 

satishes the optical selection rule. However, the same phase between A^ and B^ sublattices 
competes with the opposite phases of the and B^ ones, which predominantly reduces 
the peak intensity. Concerning the (n^ —>■ n^) and (rig —>■ n^) excitation channels, the 
= 01, = 10) and (ugUg = 01, UgUg = 10) threshold peaks are, respectively, 

revealed at {uj = 0.502 eV; oj = 0.51 eV ) and {u = 1.205 eV; uj = 1.209 eV). While the 
intensities of the latter two are comparable with the other intragroup excitation peaks, 
those of the two formers are predominant. This is because all the subenvelope functions 
of the 712 = 0 (1) initial and = 1 (0) hnal states on the A^ and B^ sublattices have 
the same (opposite) phases. In general, three categories of intragroup inter-LL excitations 
exhibit similar characteristics, including peak intensity, twin-peak structure, and selection 
rule. 

The intergroup inter-LL excitations exhibit single and twin peaks, depending on the 
quantum modes and the asymmetry of LL energy spectra. The absorption categories 
related to the second group of LLs, the —>■ n^, —)■ —)■ n^, and —)■ 

transitions, present solely single-peak structures. For the two formers, the initial and hnal 
states of available excitations satisfy e.g., 00, 00, 11, 11, and others (Fig. 

4(a)). Especially, their absorption spectra begin, respectively, with the peaks of n\n 2 = 11 
{u = 0.25 eV; marked by green arrow) and = 11 (cu = 0.245 eV), according to the 
unusual ordering of quantum numbers in the second group (Fig. 3). And then the threshold 
peaks, n\n 2 = 00 and nlnl = 00, are merged together at w = 0.265 eV and become a 
higher-intensity single peak. Concerning the ?t ,2 —>■ Tig and rig —)• excitations, the inter-LL 

transitions with — = 2 and 0 can induce the single-peak structure. The former are 
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much higher than the latter, mainly owing to the special relations among the subenvelope 
functions. For instance, the 20 single peak satishes the optical selection rule where the 
quantum modes in the A^, A^, sublattices of = 2 are, respectively, identical to those 
in the and A^ sublattices of Ug = 0. Furthermore, the corresponding sublattices 

with the same phases are responsible for the predominant peak intensity, compared to the 
weak peaks of 00, 11, and so on. On the other hand, the —)■ nl and —>■ Ug excitations 
exhibit the twin-peak structures. This is due to the fact that the inter-LL transitions obey 
the selection rule of = 1, similar to that of the intragroup ones. The twin peaks 

possess very different intensities in which one resembles the intragroup peaks while another 
is comparable to the intergroup ones. For example, the 21 and 12 peaks are, respectively, 
weak and strong as indicated by red arrows in Fig. 4(b). Exceptionally, the threshold 
peak 01 is a single one since the 10 transition is forbidden by the Fermi-Dirac distribution 
function. 

The splitting of LLs results in critical differences between magneto-optical absorption 
spectra corresponding to 4/6 (Figs. 5(a)-5(e)) and 2/6 localization centers in terms of in¬ 
tensity, frequency, and structure of peaks. This is mainly owing to the abnormal quantum 
number ordering and the signihcant side modes during the LL anti-crossings [^. Ac¬ 
cordingly, the dissimilarities of two spectra are clearly revealed in the inter-LL excitations 
related to the low-lying LLs. For the threshold peak 01, that of the former has much higher 
intensity and frequency {uj ~ 0.15 eV) (Fig. 5(a)). There are a lot of low-frequency single 
peaks when the magnetic held is within the range of LL anti-crossings. For example, the 
= 0 LL also has the non-negligible side mode of three zero points, being responsible 
for the anti-crossing of the = 0 and = 3 LLs and the appearance of the single 
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peaks 02 and 04. Furthermore, many pairs of LLs in the anti-crossings such as = 1, 
_ 2 ^ gQ present several extra single peaks. Conse¬ 

quently, similar characteristics also occur in the initial absorption peaks of the intergroup 
transitions between the hrst LL group and others. Different from the spectrum at 2/6 
center, the —>■ rig and rig —>■ excitations only present the inter-LL transitions with 

\n 2 ^ — ng’'^! = 0. In short, the LL anti-crossings play an important role in distinguishing 
the magneto-absorption spectra due to the splitting LLs. 

Specially, the magneto-absorption spectrum of ABA-stacked trilayer graphene is worthy 
of a detailed consideration for the wide-frequency range. Since LLs are four-fold degenerate, 
the inter-LL excitations at 4 different localization centers, 1/6, 2/6, 4/6, and 5/6, are 
identical. Only hve categories of absorption peaks are revealed in the absorption spectrum. 


as shown in Figs. 6(a)-6(c). They correspond to those of AB-stacked bilayer 43,44 and 


monolayer graphenes 43 , mainly owing to the absence of the inter-LL transitions related 


to a pair of linear bands and two pairs of parabolic ones. Two categories of intergroup 
excitations, 2*^ —)• 3^ and 3*^ —>■ 2'^, exhibit twin peaks due to the asymmetric LL energy 
spectrum. The intragroup inter-LL excitations of the second and third LL groups are 
mostly present in twin-peak structures except for the threshold 01 peak. The absence of 
10 transition is due to the fact that the n = 0 LL is occupied. Especially, the category 
of 1*^ transitions generally displays single peaks with dominant uniform intensity 

because of the symmetric LL energy spectrum of the quantized linear bands. There is only 
a twin-peak structure for 21 and 12 (red arrows in Fig. 6(b)), reflecting the distortion of 
the linear bands near the Dirac points. Moreover, the peak frequencies are proportional to 
■y/^, similar to those in three absorption categories of the AAA stacking. On the other 
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hand, there is no simple relationship in the other inter-LL excitations of the ABA, ABC 


and AAB stackings. 


The stacking conhgurations can diversify the magnetic quantization and thus the magneto¬ 
absorption spectra. Four typical trilayer graphenes, AAB, AAA, ABC, and ABA stack¬ 
ings, are quite different from one another in the excitation category, structure, intensity, 
and number of absorption peaks. The differences mainly come from the special relations 
among the subenvelope functions on the six sublattices. The AAA trilayer graphene, being 
regarded as three separated monolayer graphenes, only has three categories of intragroup 


absorption peaks 38,39 . The main reason for the absence of the intergroup excitations 


is that all the LL wavefunctions are the linear symmetric or asymmetric superpositions of 
the subenvelope functions in different layers. The absorption spectrum only exhibits the 
single-peak structures with uniform intensity except for the threshold peak in each category 
with half of intensity. On the other hand, the ABC-stacked trilayer graphene possesses nine 


categories of absorption peaks with nonuniform intensities 40,41 , similar to those of the 


AAB stacking. The whole spectrum, excluding the threshold peaks, exhibits the twin-peak 
structure. Overall, the intragroup excitations are stronger than the intergroup ones. This 
can also be observed in the ABA stacking (Fig. 6(a)-6(c)), being in sharp contrast with the 
AAB stacking. Concerning the number of absorption peaks, it is higher for the ABC and 
AAB stackings but lowest for the AAA one. The critical differences among the four stacked 
trilayer graphenes can provide the useful information in identifying the specihc stacking 


configurations by the experimental measurements 54-59 . 


4. Conclusion 
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The magneto-optical properties of the AAB-stacked trilayer graphene have been investi¬ 
gated using the generalized tight-binding model. This method can also be extended to study 
the other essential physical properties of few-layer graphenes and other 2D materials. An 
optical selection rule that the sublattice of the initial state has the same quantum 

mode with the (A^) sublattice of the hnal state is deduced from the spatial distributions 
of LL wavefunctions. The main features of magneto-absorption spectra, including the ex¬ 
citation category, intensity, structure, frequency, and number of absorption peaks, are fully 
explored in a wide frequency range. There are nine categories of inter-LL optical transitions 
in which the intragroup and intergroup ones quite differ from one another. The magneto¬ 
absorption spectrum could be enriched by the splitting LLs, arising from the destruction of 
inversion symmetry in the AAB stacking. Both single- and twin-peak structures are clearly 
revealed in the magneto-absorption spectra. The main characteristics of LL wavefunctions 
and energy spectra and the Fermi-Dirac function are responsible for the above-mentioned 
characteristics. The predicted feature-rich magneto-absorption spectra with nine excitation 


categories can be further verihed by optical spectroscopy techniques 54-59 . 


The intergroup and intragroup absorption peaks are signihcantly different in intensity, 
structure, and frequency. In general, the former are much higher than the latter, depending 
on the phase and amplitude relations between the subenvelope functions on six sublattices. 
Due to the quantum modes and asymmetric LLs, the former have single and twin peaks 
while the latter display twin-peak structures except for threshold ones. The energy spacings 
between the twin-peak structures are much wider for the intergroup excitations compared 
with the intragroup ones. The critical differences in spectra at 2/6 (5/6) and 4/6 (1/6) 
localization centers mostly he in the excitations associated to the low-lying LLs. This is 
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caused by the significant side modes during the LL anti-crossings. A lot of low-frequency 
extra single peaks are present in the latter when the magnetic field is within the range 
of LL anti-crossings. The inter-LL excitations of the higher-lying LLs are quite similar in 
certain respects for different localization centers. 

The AAB-stacked trilayer graphene exhibits the unique magneto-optical properties 
among four typical stacking systems. It possesses the most plentiful magneto-absorption 
spectrum, mainly owing to the splitting LLs and the frequent LL anti-crossings. The AAA 
stacking, the superposition of three monolayer graphenes, has three categories of intragroup 
excitation peaks. The magneto-absorption spectrum only has many single peaks with a 
uniform intensity except for the threshold ones. The ABA stacking, being regarded as the 
combination of monolayer and AB-stacked bilayer graphene, only presents five categories of 
absorption peaks. Only peaks arising from the slightly distorted Dirac cones show the uni¬ 
form spectral intensity in the single-peak structures. Two pairs of parabolic bands induce 
two intragroup and two intergroup excitations categories, in which the former are stronger 
than the latter. All the peaks, the threshold ones excepted, have the twin-peak structures. 
The ABC stacking possesses the same number of intra- and inter-group absorption peak 
categories as the AAB one. However, the intensity relation between the intragroup and 
intergroup peaks is similar to that of the ABA stacking, but in sharp contrast with the 
AAB one. 
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FIGURE CAPTIONS 


Fig. 1 - The interlayer atomic interactions and the geometric structure under a uniform 
magnetic held Bqz (a). The shaded region corresponds to a rectangular unit cell. The hrst 
and second layers have the same (x, y) projections, (b) The low-energy band structures 
of AAB-stacked trilayer graphene with (c) a narrow energy gap. The solid and dashed 
curves represent the results by the use of the tight-binding model and the hrst-principle 
calculation, respectively. 

Fig. 2 - The distributions of the low-lying LL wavefunctions of AAB-stacked trilayer 
graphene with the six distinct sublattices centered at the 2/6 localization under Bq = 40 
T. The unit of the x-axis is m/2RB, where m represents the m-th A or B atom in the 
enlarged unit cell. 

Fig. 3 - The distributions of the conduction (a) and valence (b) LL wavefunctions 
of AAB-stacked trilayer graphene with the six distinct sublattices centered at the 2/6 
localization under Bq = 40 T. The unit of the x-axis is m/2RB, where m represents the 
m-th A or B atom in the enlarged unit cell. 

Fig. 4 - Magneto-absorption spectra of AAB-stacked trilayer graphene in (a), (b), (c), 
(d) and (e) near 2/6 center with respect to different frequency ranges at Bq = 40 T. The 
numbers with and without the underline stand for the quantum numbers of the occupied 
and unoccupied LLs, respectively. The blue, green, and red colors correspond to the first, 
second, and third groups of LLs. 

Fig. 5 - Magneto-absorption spectra of AAB-stacked trilayer graphene in (a), (b), (c), 
(d) and (e) near 4/6 center with respect to different frequency ranges at Bq = 40 T. The 
numbers with and without the underline stand for the quantum numbers of the occupied 
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and unoccupied LLs, respectively. The blue, green, and red colors correspond to the hrst, 
second, and third groups of LLs. 

Fig. 6 - Magneto-absorption spectra of AAB-stacked trilayer graphene in (a), (b), 
and (c) near 2/6 center with respect to different frequency ranges at Bq = 40 T. The 
numbers with and without the underline stand for the quantum numbers of the occupied 
and unoccupied LLs, respectively. The blue, green, and red colors correspond to the hrst, 
second, and third groups of LLs. 



M <- K (iM)-> r k^{2n/^b) M <- K (1/4)-> r 


Figure 1: The interlayer atomic interactions and the geometric structure under a uniform 
magnetic held Bqz (a). The shaded region corresponds to a rectangular unit cell. The hrst 
and second layers have the same (x, y) projections, (b) The low-energy band structures 
of AAB-stacked trilayer graphene with (c) a narrow energy gap. The solid and dashed 
curves represent the results by the use of the tight-binding model and the hrst-principle 
calculation, respectively. 
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Figure 2: The distributions of the lo'w-lying LL 'wavefunctions of AAB-stacked trilayer 
graphene with the six distinct sublattices centered at the 2/6 localization under Bo = 40 
T. The unit of the x-axis is m/2RB, where m represents the m-th A or B atom in the 
enlarged unit cell. 
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Figure 3: The distributions of the conduction (a) and valence (b) LL wavefunctions of AAB- 
stacked trilayer graphene with the six distinct sublattices centered at the 2/6 localization 
under Bq = 40 T. The unit of the x-axis is m/2RB, where m represents the m-th A or B 
atom in the enlarged unit cell. 
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Figure 4: Magneto-absorption spectra of AAB-stacked trilayer graphene in (a), (b), (c), 
(d) and (e) near 2/6 center with respect to different frequency ranges at Bq = 40 T. The 
numbers with and without the underline stand for the quantum numbers of the occupied 
and unoccupied LLs, respectively. The blue, green, and red colors correspond to the first, 
second, and third groups of LLs. 
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Figure 5: Magneto-absorption spectra of AAB-stacked trilayer graphene in (a), (b), (c), 
(d) and (e) near 4/6 center with respect to different frequency ranges at Bq = 40 T. The 
numbers with and without the underline stand for the quantum numbers of the occupied 
and unoccupied LLs, respectively. The blue, green, and red colors correspond to the first, 
second, and third groups of LLs. 
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Figure 6: Magneto-absorption spectra of AAB-stacked trilayer graphene in (a), (b), and (c) 
near 2/6 center with respect to different frequency ranges at Bq = 40 T. The numbers with 
and without the underline stand for the quantum numbers of the occupied and unoccupied 
LLs, respectively. The blue, green, and red colors correspond to the hrst, second, and third 
groups of LLs. 
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